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Abstract 



We study the dynamics of electroweak bubbles in the scenario of a strong first order inverse 
u ; 

electroweak phase transition at the TeV scale involving the global structure of the nonlinear sigma 



field in the littlest Higgs model with T parity. Employing the onedoop order finite temperature 
effective potential, we find that the pressure in the symmetric phase i.e., inside the bubble is 



— always greater than that in the asymmetric phase i.e., outside the bubble, so that the bubbles are 

qJ expanding. By calculating the fluid velocities in the two phases we arrive at the condition of a 

D ■ 

supersonic deflagrated motion of the bubble walls. We then discuss the generation of gravity waves 



from the collisions of such bubbles as well as from the turbulence of the plasma. 
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I. INTRODUCTION 

In recent years, dynamics of electroweak bubbles, associated with a first order phase 
transition, has been studied in the Standard Model (SM) [1-11] and its extensions [12-14], 
Two-Higgs Doublet Model (TDHM) [15], the Minimally Supersymmetric Standard Model 
(MSSM) [16, 17] as well as in model-independent way [18,19]. Also, aspects of gravity waves 
(GW), as would be generated by bubble collisions and due to turbulence in the plasma have 
been investigated [20-26]. 

Knowledge of bubble dynamics and the concomitant CP violation [27] help us to un- 
derstand the electroweak phase transition (EWPT) and electroweak baryogenesis (EWBG) 
scenario in more detail in a specific model. 

Depending on the strength of the first-order phase transition and the nature of the finite- 
temperature effective potential (FTEP) in a model calculation, the motion of the bubble wall 
may belong to various categories, viz., deflagrations, detonations, hybrids or runaway [19]. 
Also, the motion can be supersonic, Jouguet or subsonic depending on whether the velocity 
of the plasma inside the bubble is greater than, equal to or less than the sound velocity in the 
medium, respectively. Usually, deflagrations are subsonic and detonations are supersonic. 
However, no clear-cut classification can be made in this regard. Although hydrodynamical 
equations give general features of bubble wall motion, microscopic analysis including friction 
of the bubble wall is necessary for getting the detailed information regarding the bubble wall 
velocity [13, 19]. 

Gravity waves, which are considered to be the essential features of Einstein's General 
Theory of Relativity, have not been experimentally detected so far, although efforts are on 
in this direction [26]. Gravity waves at the electroweak scale are believed to have been 
generated either by bubble wall collisions or by turbulence of the plasma. 

The littlest Higgs model (L 2 HM) [28] and the littlest Higgs model with T parity (LHT)[29- 
31] are economical beyond Standard Models (BSM) which can solve the little hierarchy 
problem. The LHT conforms to the electroweak precision data. In order to explore the 
prospects of EWPT and EWBG, finite-temperature calculations have been performed [32- 
34] in these models. In Ref. 34, an inverse strong first order EWPT has been observed in 
the global structure of the effective potential at the value of the physical Higgs field, h — 1.1 
TeV with transition temperature, T c = 0.9 TeV. 



In Ref. 35, a two-step baryogenesis scenario associated with the inverse phase transition 
in the TeV scale has been presented. In the first step, the Universe makes a transition from 
an electroweak broken phase above T c to a symmetric phase below T c . Bubbles of symmetric 
phase are formed in the background of asymmetric phase. Baryon number violations take 
place within these bubbles due to sphaleron transitions induced by T even massless gauge 
boson fields. In the second step, a cross-over takes place at T ps 0.1 TeV and the process of 
baryogenesis gets completed between the temperatures, 0.9 TeV and 0.1 TeV. 

In view of the proposed new aspects [34] of the EWPT at the TeV scale and the associated 
EWBG scenario [35] , we expect to find new features of the bubble dynamics as well as GW 
in the LHT. 

The purpose of this paper is to examine the properties of the symmetric phase bub- 
bles associated with an inverse strong first order EWPT and a new baryogenesis scenario 
therewith. In section II, we demonstrate the expansion of the symmetric phase bubbles. In 
section III, we calculate the bubble nucleation rate and time and in section IV, we study 
the bubble wall velocity. In section V we investigate the features of gravity wave generation 
in our model. Finally, in section VI we make some concluding remarks. 

II. EXPANSION OF THE SYMMETRIC PHASE BUBBLES 

The presence of bubbles in the model under consideration is schematically shown in the 
Fig.l. Here we have a situation which is opposite to that in the case of SM, viz., we have 
the symmetric phase bubbles in the background of the asymmetric phase. Whether these 
bubbles will contract or expand will depend on the difference of pressure in the two phases. 
The pressures are completely determined by the FTEP [34] and can be written as [4], 

Pin = -Veff (h = 0,T) , Pout = -V eff (h = l.lTeV,T) , (1) 

where, pi n and p out are the pressures inside and outside the bubbles respectively. It should 
be noted that h = 1.1 TeV is the value of the physical Higgs field, where the phase transition 
is observed. 

So far as the FTEP is concerned, two comments are in order here. Firstly, in the cal- 
culation of FTEP in Ref. 34, the effects of mirror fermions, which are instrumental in the 
non-minimal CP violation [36] in the LHT models, were not taken. Subsequently, we have 
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FIG. 1. Electroweak bubbles in the inverse phase transition scenario in the littlest Higgs model 
with T parity. Inside the bubbles, the VEV of the physical Higgs field, < h >= 0. Outside the 
bubbles, < h >= 1.1 TeV. 



found that the inclusion of these fermions in the FTEP does not change the features of the 
phase transition - the value of T c changes from 0.925 TeV to 0.905 TeV and the position of 
the minimum of the Higgs field remains the same. However, to be accurate, for the pres- 
sure calculations here, we have included the mirror fermions in the FTEP. The value of the 
mirror fermion coupling constant k has been taken to be 0.7 consistent with a reasonable 
value of the mirror fermion mass, viz., 0.5 TeV [37]. The second comment concerns the 
recently discovered [38, 39] SM Higgs of mass, m# ~ 125 GeV at the LHC. This translates 
into value of Higgs quartic coupling constant, A = [jriH ^"> = 0.13, v being the SM Higgs 
VEV, 246 GeV. With this value of A, we get a set of UV completion factors as a = —0.01, 
a' = —0.0002 for our FTEP, consistent with the experimental value of the SM Higgs VEV 
at zero temperature. It may be noted here that the general procedure for obtaining the UV 
completion factors in our model is described in Ref.34. 

With the above considerations and Eq.l we have calculated the quantity, (p in — p out ) of 
which the plot is shown in Fig. 2. We observe that the excess pressure within the symmetric 
phase bubbles is always positive. From this we may conclude that after the inverse EWPT at 
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FIG. 2. Excess pressure within the bubbles of the symmetric phase as a function of temperature. 

T c =0.9 TeV, the symmetric phase bubbles will expand in the background of the asymmetric 
phase. 

III. NUCLEATION RATE AND TEMPERATURE 



The excess free energy in a true vacuum bubble, which is a symmetric phase bubble in 
the present case, can be written as, 



AF(T) =-- 4vr / i 
'o 



KS'i +^ T) 



dr. 



(2) 



where, we have assumed that the bubble is spherical and R is the radius of the bubble. 
The derivative term in Eq.2 corresponds to a surface energy and the potential term a vol- 
ume energy [3]. In the thin wall approximation, appropriate for a strong first-order phase 



transition, the derivative of the Higgs field can be written as [40], ±^2V (h,T) in the limit 
of exact degeneracy of the potential minima at T = T c . A positive value of the derivative 
would reflect the fact that in our case as we go from the centre of a symmetric phase bubble 
to the boundary, the VEV of h increases from to 1.1 TeV. To avoid using an ansatz of the 
field in the r space, we then change the integration variable from r to h and approximate 
the surface term as, 47r_R 2 cr(T), where cr(T) = J h ~ Q ' ^/2V (h,T)dh can be called the surface 
energy density. We also write the volume term in Eq.2 as ^nE?V (h, T) where V (h, T) is 
an average effective potential inside the bubble or minus of an average pressure inside the 



bubble which can be approximated as the minus of the excess pressure inside the bubble, 
since Eq.2 defines an excess quantity. Thus, we write the volume term as, — ^irR 3 Ap(T). 

Now, we get the radius of a critical bubble by minimizing AT (T) from the equation, 
with AT(T) = 4xR 2 <t(T) - ^R 3 Ap(T) and thus R C (T) 

The bubble nucleation rate per unit time per unit volume can be written in terms of the 
excess free energy ATc(T) of the critical bubble as [41,17], 
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(3) 



The bubble nucleation temperature T/v, which is somewhat lower than T c in the case of 
first order phase transition, is defined as the temperature at which the rate of nucleation 
of a critical bubble within a horizon volume is equal to the Hubble parameter at that 
temperature. Since the horizon scale is approximately /J(T)" 1 the nucleation temperature 
can be defined by the equation, 



H(T N ) ~ 1.66 V&(T; 



■N 



rp2 

2 N 



H d (T N ) mpi 

where g*{T) is the relativistic degrees of freedom and mpi ~ (1.22 x lO 19 ^^^) is the Planck 
mass. From Eqs. (3) and (4), we get, 
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(5) 



Eq.5 can be graphically solved for T/v- Denoting the LHS of Eq.5 as /(Tjy), we show in 
Fig.3 a plot of /(T/v) vs. T/v. The graph shows that /(Tv) ~ for T/v ~ 0.5 — 0.6 TeV. 
The nucleation temperature can also be determined from the consideration of the fact that 
the space-time integrated nucleation rate (i.e., the nucleation probability) should be of the 
order of 1 [13]: 



V N {t)V c dt 



(6) 



tc 



Now the horizon scale is roughly 2t and so, volume = (2t) 3 . Using the time-temperature 
relations [40, 13], 



t = 0.301g; 



-l/2 m Pl 
y2 ' 

— = -HT 
dt 



(7) 



we get the integral in Eq.6, which is the total nucleation probability of a single bubble, as, 
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FIG. 3. A plot of /(T/v) vs. T/v where /(Tv) is the function on the left hand side of Eq.5 and T/v 
is the bubble nucleation temperature. 

Results of evaluation of this integral for different values of T/v are shown in Table 1 . We see 
from Table 1 that consistent with Eq.6, the value of T/v should be 0.635 TeV. 

IV. THE NATURE OF BUBBLE AND WALL VELOCITY 



The entropy, energy and enthalpy densities are defined in terms of the pressure p(T) as, 



s(T) = (dp(T))/dT = p'(T), 

p(T)=Ts(T)-p(T), 

u(T) = p(T)+p(T)=TS(T). 

We have then the following equations, valid for a first-order phase transition. 

P+(T C ) =P-(T C ), 

S+ (T c )^s_(T c ), 

P+ (T c )^p_(T c ), 

and ij + (T c ) ^ U- (T c ) , 



(9) 



(10) 



where, the quantities with the subscript, '+' refer to those in the high-temperature phase 
( i.e., the phase outside the bubble) and the subscript, '— ' is for the quantities in the low- 
temperature phase (i.e., the phase inside the bubble). 



TABLE I. Nucleation temperature, excess free-energy and the nucleation probability (integral I in 

Eq.8). 
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The latent heat L is defined as the difference between the energy densities of the two 
phases at T = T c and thus, 



L = T c \p'(T c )-p'_(T c )} 



(11) 



Solving hydrodynamical equations [19,42,43], we get, 



w+7+ v 1 + p+ = u^ 2 _v 2 _ + p_ 



(12) 



where, 7 = 1/vl — v 2 , v being the velocity being assumed to be aligned in a particular 
direction, say, the z-direction. 

Equivalent ly, we get, 



v + v. 



P+-P- ^± = P- + P+ 
P+-P- v - P+ ~ P- 



(13) 



whence, we obtain, 



'(p+- 


-p-)(p- 


+ P+) 


(p+- 


- p-)(p+ 


-p-Y 


l(p + 


-P-)(P^ 


--P-) 



(14) 
(p+-p-)(p-+p+) 

where, it is implied that, p + = p+(T + ), p_ = p_(T_) etc., T + being the temperature outside 
the bubble and T_ inside. 

The bubble motions can be broadly classified as [13], 

(i) Deflagration : The velocity of the plasma inside the bubble is more than outside, 
v- > v+ , 

(ii) Detonation : The velocity outside is more than inside, v + > i>_. 

Deflagration may further be subdivided as, strong: v _ > c s , Jouguet: V- = c s and weak: 
v _ < c s . Here, c s is the velocity of sound inside the bubble, which has the value l/\/3 
for a relativistic plasma [13]. The strong deflagration is therefore supersonic and the weak 
deflagration is subsonic. A detonation will be called strong if v _ < c s , Jouguet if i>_ = c s 
and weak if v _ > c s . 

Deflagrations and detonations are also characterized by the following temperature con- 
straints [13, 18, 19]: 

Deflagration : T + > T N > T_, 
Detonation : T_ > T + = T N . (15) 

The above characterizations show that the detonation cases are easier to investigate because 
of the condition, T + = T N . 

Now, coming to the case of LHT, let us first assume that the bubble motion is detonated, 
so that T + = Tjy=0.635 TeV. Fig.4 shows the result of calculation of v + and V- as a function 
of T_ for the parameter space / = 500 GeV, ttih =125 GeV, k = 0.7 and T + = T/v (= 0.635 
TeV). The plot clearly shows that v_ > c s > v+, which is in contradiction with the velocity 
characterization of detonation. Hence we can rule out the case of detonation in the case of 
LHT. Also as i>_ > c s always, we are getting a case of supersonic deflagration [18] for the 
bubble wall motion here. 

To study the bubble wall motion, let us use the following so-called bag equations of state 
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FIG. 4. Assumming detonation motion (T+ = Xjy ), plot of velocity of plasma outside the bubble: 
v + (lower plot: blue-dashed), sound velocity: cs(= -?=), middle green straight line) and velocity of 
plasma inside the bubble: u_( upper plot: blue-dotted). 



[13]: 



p+ = a + T 4 + s, p + = -a + T 4 — e, p_ = a^T A , p. 



-a-T 4 



and the entropy density, 
where, 



S± = o°± T3 ' 



a± = ^9*(T±) 



(16) 



(17) 



(18) 



are numbers related to the number of relativistic species present in the plasma. It is inter- 
esting to see the physical significance of the quantity e. Since at T c , the pressure difference 
between the two phases is zero, we get from Eq.16, 



a-)T? 



3e 



Then, from the expressions of entropy and latent heat we get, 



L = Ae 



(19) 



(20) 



which shows that e is closely related to the amount of latent heat L of phase transition. 
We now define two important quantities upon which the bubble wall velocity depends: 

e e 



a 



a + T 



■ l • 
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a + T N 



(21) 
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From hydrodynamics, we get a relationship between the plasma velocities in terms of a as 

[44], 



6v- ~ l ~ 2 V V 6u_ "^ 2 



*•'+ 



a z 



!« 



l + a 



(22) 



where, the +(— ) sign corresponds to the detonated (deflagrated) motion of the bubble. 

As we have argued earlier, the motion of the bubble is deflagrated rather than detonated. 

Accordingly, we shall study, henceforth, the deflagrated motion of the bubble, taking the 

'— ' sign in Eq.22. 

A deflagrated bubble proceeds via the generation of shock waves. The expressions of the 

velocities v\ and v 2 of the fluid before and behind respectively of the shock front in the frame 

of the shock front are [13], 

1 /3Tf+T^\ 1 

In the laboratory frame the fluid within the bubble, i.e., behind the bubble wall is at rest, 
as the bubble is isotropically expanding, and also ahead of the shock front, where the 
temperature is Tjy, which belongs neither to the low-temperature nor to the high temperature 
phase, the fluid is at rest. Thus, the wall velocity in the laboratory frame, 

v w = \v-\ (24) 

because, i>_ is the velocity of the fluid inside the bubble in the reference frame of the wall. 
Thus in terms of the bubble wall velocity also, it can be said that the motion is supersonic 
if we can show that V- > c s for all temperatures T + outside the bubble. The fluid velocity 
Vf between the two fronts can be written either in terms of Vi, v 2 or v + , t>_: 

Vi -v 2 v+ — V- 



1 — f if 2 1 — V + V- 

i/3(ajv — a) 



v .f 

Putting the values of V\ and v 2 we get, 

v + — V- 

1 - v+v- a/3 (a N + a) (3a + a N ) 
Eliminating v + from Eqs.25 and 26 we get an expression of V- in terms of a and a^'- 

i>_ = (V3(a + — aAr)(3o; + + a7v)(a + + 3«Af)) _1 

3«Ar)(a+ — 2a+(l + 4a + )aN + &%) 



(25) 



(26) 



a/(3« + + «Ar)(a; 

J r2^f((?>a + + aAr)(o; + + 3cxn) (8a+a% + a% — 2a 
y +a\{l + 4a JV (-l + 4a N ))) 



+ a N 
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'l + 2a 



N 



(27) 
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In LHT, g*{T > 500 GeV) ~ 214, 7^=0.635 TeV, T c =0.91 TeV, L=0.2464 TeV. Then we 
get, a N = 0.00538.. 

Keeping the latent heat and therefore a^ constant, we calculate w_ varying a. A plot of 
v _ against a (Fig. 5) shows that the motion of the fluid inside the bubble, and therefore the 
bubble wall velocity v w in the laboratory frame is supersonic. We can also plot V- against 
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FIG. 5. Plot of V- against a for fixed a^ following Eq.27. For various values of the temperature 
T + , V- is always greater than the sound velocity, c s =l/\/3=0.577. 

a for various values of v + following Eq.22 instead of Eq.27 where v + has been eliminated. 
Such graphs are shown in Fig. 6. It is seen from these graphs that for various values of v + 
also the supersonic nature of wall velocity i.e., v w = |t>_| > c s is maintained. 

V. GRAVITY WAVES 

Gravity waves [46-48] may be generated at the time of the first-order EWPT by two 
mechanisms [20-26]: (i) Collision of bubble walls, (ii) Turbulence created by the bubble 
expansion. Of these two processes, the latter is more effective in producing GW than the 
former and a strong first-order phase transition is conducive for this process. Since the strong 
first-order phase transition is not possible in the SM, GW cannot be generated within the 
framework of the SM. 

The frequency and energy density of GW produced in the early Universe change due to 
the expansion of the Universe. The frequency varies as a -1 and the energy density a -4 where 
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FIG. 6. Contour plot of V- against a for various values of u+. Contours denote the values of v+. 



a is the scale factor. Thus the red-shifted frequency we observe today is given by, 

/o = /* — ) 
a 



(28) 



where, a is the scale factor at the present time and a*(/*) is the scale-factor (frequency) at 
the early Universe, i.e. at the time of the origin of the GW. 

Now, because of the adiabatic expansion of the Universe, we can write, 

Substituting, T = 2.TZK = 2.35 x l(T 13 GeV, g = 3.91 we get, 

'100\ 1/3 flOOGeV\ 



(29) 



- = 8.0 x 10- 16 



(30) 



ao \ 9* J V ^* / 

The wavelength of GW in the early Universe should be a fraction of the Hubble size H^ 1 . 



Hence, it is instructive to write f in terms of /*/-H*. Thus, we write, 

where, H* is given by the Friedmann equation, 

tf 2 - 8nG n 



(31) 



(32) 
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p* being the radiation energy density, 



P* = -go- (33) 



From the above equations, we get a suggestive expression of frequency of GW in the present 

Universe, 

/o = 1.65 x 10- 5 f-^V /6 (^-) (— ^ Hz, (34) 

J V100/ \hJ \l00GeVJ ' V ; 

where, we have used the numerical values [40], mpi = 1.22 x 10 19 GeV and lGeV^ 1 = 
6.5822 x 10- 25 sec _1 . 

Let us first consider GW generation from bubble wall collisions. During EWPT, the 
bubble nucleation rate per unit time and unit volume may be written as [26, 49], 

r N (t) = r (t,)e^-^, (35) 

where, ti is an initial time and (3 is a parameter setting the time and length scales of the 
phase transition as /3 _1 and v m fi~ x respectively, v w being the bubble wall velocity.Assuming 
/* ~ j3, we can write from Eq.34, 

C l = 1.60 x 10" 5 (-* ) V6 ( T * ) (I-) Hz. (36) 

J0 vioo/ \iooGevj \ H *J 

A recent simulation [20, 26] gives the dependence of the peak frequency, j^k on v w as ; 

* — > * - 5 (u=^s) &T (iooM (I) -• <"> 

which is interestingly similar to Eq.36. We may also observe that the peak frequency of GW 
due to bubble collisions does not depend very sensitively on the bubble wall velocity. 

The fractional gravitational energy density (or intensity) in the early Universe is defined 

as, 

&gw* — (38) 

Ptot* 

and its relation with the corresponding value fi^vU a ^ the P resen t Universe is given by, 

Using the value of the Hubble constant in the present Universe [40], H = 2.1332/i x 
W- 42 GeV(h ~ 0.7), Hi = (8tt 3 /90) Gg*T* = 1.85 x 10- 38 #*T 4 GeV~ 2 and Eq.30, we get, 

h 2 n^ = 1.67 x 10- 5 x I J Q^- t . (40) 
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Eq.40 gives a relation between the intensities of the gravity waves in the early Universe and 
the present Universe. 

The GW energy density in the early Universe is given by [49], 

«-=* 2 fc) 2 (xK' (41) 

where, the efficiency factor (k) is defined by the equation, 

p k = kap rad , (42) 

oiprad being the latent heat and pk is the kinetic energy density. 
We have here, 

Put* = (1 + a)p rad (43) 

From Eq.40 we then get, 

tfngft. = 1.67 x 10- (f) "V ( T ^) 2 [f) \ ,44) 

A recent fit gives [20, 26], 

^,,4 1Xl 0- X (M) 1 "V(^) 2 (f) 2 (^) («, 

which is close to Eq.44. 

The GW produced by the stirring of the plasma or by the turbulent bulk motion of the 
plasma is somewhat different from that due to bubble collisions. In this case the relevant 
length scale is the so-called 'stirring scale' [26], which is approximately given by L s ss 2R b , 
where R b is the bubble radius. For the largest bubble, R b ~ v w (3~ l . A better approximation 
[26] for the radius is 

R b n?>v w \nU^\p-\ (46) 

As frequency, /* ~ Lj 1 , we get from Eq.36, 

/o ^ 8xlo -«(^) 1/6 f^^U^Hz (47) 

J0 V100/ \100GeVJ \H*R b ) V ; 

For the case of GW produced by turbulence, the expression of the peak intensity of the wave 
can be cast in the form [26, 50-53], 



100\ 1/3 / a \ 3/2 / (R b Hl 
gj \ (l + a)J ll + 4 7 3 -' 



^ 2 «^3.5xl0- 5 ( — ) [kj^] I n y r°"H 1 (48) 
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We now turn to the calculations of the frequency and intensity of GW in our model. We 
have written the bubble nucleation rate both as functions of temperature (Eq.3) and time 
(Eq.35). The rates are equal at time and temperature related by Eqs.7. From Eq.35, we 
get, 



fi = ^ (49) 

J- TV 



Then, from Eqs.7 and 3 we obtain, 



| = T-L (AF C (T)/T) (50) 

where, we have neglected the time variation of the prefactor in comparison to that of the 
exponent as the nucleation rate is mainly driven by the exponential factor [26]. 

To get the derivative T-^ (AFc(T)/T) at various temperatures in the desired range we 
have made a polynomial fit of AFc(T)/T up to fifth power of T from its values calculated 
at a number of temperatures and tabulated in Table 2. The fitted polynomial is, 

9\J. = -3877.96+87751.84T-572424.06T 2 +1584342.13T 3 -1952687.4T 4 +881672.71T 5 

(51) 
Differentiating the function (Eq.51) with respect to T we can get the value of /3/H at a 
particular temperature. 

In Ref.19 an expression of efficiency factor for supersonic deflagration has been derived 
in model-independent way: 

k = fa + (yw - c 3 ) x 5 + ( — — — j (fa - fa - (cj - cs) x 5), (52) 

\cj — Cs / 

a 2/5 , Jo. , / Jol 



fa = ^7F, fa = v - 5 = -0.9 In v ^ (53) 

0.017+ (0.997 + af 5 0.135 + V0.98 + a \l + v^, 

and cs and cj are the sound and Jouguet velocities respectively. 



CJ = ^vs+v^+v 3 (54) 

1 + a 

In Fig. 7, we show the variation of GW signal generated by bubble wall collisions against 
temperature as well as the bubble wall velocity. The plot shows that the intensity is prac- 
tically independent of the bubble wall velocity. This behavior can be understood from the 
v w dependence of the intensity, as demonstrated by Eq.45. This feature of velocity inde- 
pendence has been found here in the case of frequency also as we have already mentioned. 
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TABLE II. Values of various quantities, which are useful for studying bubble dynamics. 



T(TeV) 


a (TeV 3 ) 


R c (TeV- 1 ) 


AF c (TeV) 


AF C /T 


0.10 


0.06198 


17.48 


79.28 


792.87 


0.15 


0.06197 


17.46 


79.15 


527.68 


0.20 


0.06196 


17.41 


78.70 


393.50 


0.25 


0.06187 


17.44 


78.90 


315.62 


0.30 


0.06170 


17.66 


80.63 


268.77 


0.35 


0.06153 


17.89 


82.54 


235.85 


0.40 


0.06130 


18.19 


85.02 


212.56 


0.45 


0.06102 


18.60 


88.48 


196.63 


0.50 


0.06067 


19.19 


93.59 


187.19 


0.55 


0.06023 


20.00 


100.91 


183.48 


0.60 


0.05964 


21.12 


111.53 


185.88 


0.65 


0.05883 


22.84 


128.63 


197.89 


0.70 


0.05779 


25.41 


156.32 


223.31 


0.75 


0.05634 


29.53 


205.89 


274.52 


0.80 


0.05435 


37.14 


314.21 


392.76 


0.85 


0.05145 


56.13 


679.04 


798.87 


0.90 


0.04667 


172.06 


5786.28 


6429.20 



Next, we have compared the GW intensities for the cases of bubble collisions and turbulence 
for a fixed value of the wall velocity (t> w =0.6) and for temperatures above the nucleation 
temperature T/v=0.635 TeV in Fig.8. We observe the general trend of the signal from tur- 
bulence to be more enhanced than that from bubble collisions. Specifically, the value of the 
peak intensity arising out of turbulence is about 10 -17 whereas that due to bubble collisions 
is 10~ 20 . The calculation shows that the turbulent motion of the plasma was more instru- 
mental in producing GW in the early Universe than the bubble collisions. We should also 
notice here the kinks of the graphs at around T=0.74 TeV which are due to sudden change 
in the slopes of the graph of excess free energy (See Fig. 3) at this temperature. In Fig.9 we 
have shown the frequency distribution against temperature for GW from turbulence as well 
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FIG. 7. Plot of variation of intensity of GW from bubble wall collisions against wall velocity v w 
and temperature T. 
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FIG. 8. Plot of Intensity of GW due to turbulence (solid line) and from wall collisions (dotted 
line). 

as bubble collisions. The peak frequencies around the region of peak intensities are only 
fractions of a Hertz. 

The weak GW signals that are obtained in our calculations are difficult to detect in the 
ongoing GW detectors including LIGO. However, there are good possibilities of detection 
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in Ultimate DECIGO [26] and BBO Correlated [46], as the sensitivities of these detectors 
conform to the frequency and intensity range obtained in the present calculation. 
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FIG. 9. Plot of frequency of GW due to bubble collisions (dotted line) and from turbulence (solid 
line). 

VI. CONCLUSIONS 

In conclusion, we have studied in the present paper some aspects of the dynamics of the 
bubbles associated with first-order electroweak phase transition within the framework of the 
littlest Higgs model with T parity. We had earlier, in this model noticed a new region of 
strong first-order phase transition and proposed a baryogenesis scenario at the TeV scale. 
Our general observation in the present work is that the bubble wall motions are supersonic 
deflagrated. We have also considered the generation of gravity waves due to bubble collisions 
as well as turbulence in the plasma. As expected, the calculated intensities for the latter case 
are four order of magnitude larger than for the former case, at their peak values. Although 
the obtained intensities are quite small in magnitude, we expect them to be detected by 
some of the future gravity wave detectors. 
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